Resistivity effects in surface superconductivity of thin films in strong magnetic fields 
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Phase slips creation in the thin film in perpendicular magnetic filed with edge superconductivity is 
studied. These centers are due to thermal activation of the order parameter below superconducting 
temperature transition leading to the suppression of the superconductivity. The corresponding 
resistance is calculated. The Alsamazov- Larkin correction to the conductivity above the critical 
magnetic field destroying the surface superconductivity is studied. Such structures could be applied 
as a new system for the study of the phase slip phenomenon in one- dimensional superconducting 
wires. 

PACS numbers: 74.25.Fy, 74.25. Ha, 74.40.+k, 74.25. Op 

As it was first shown by Saint- James and de Gcnncs 
[H, superconductivity can nucleate in the thin surface 
superconducting sheath to magnetic fields higher than 
the bulk critical field H c2 < H < H c3 « 1.69iT c2 . 

Thin films reveal the most simple picture of surface su- 
perconductivity. In particular, this case was studied ex- 
perimentally in papers [2, 3], where the temperature de- 
pendencies of resistivity of thin Nb films were measured. 
Effects of the surface inhomogeneities and the sample 
shape, properties of mesoscopic size superconductors in 
surface superconductivity regime were investigated both 
theoretically and experimentally in papers 

However, it is of definite interest to consider the limit 
of the extremely thin superconducting film, where super- 
conductivity persists in the quasi- one- dimensional edge 
layer. 

It is well known that the fluctuations of the order pa- 
rameter play an important role in the physics of low- 
dimensional superconductors (thin films, wires). At tem- 
peratures above the critical temperature of the supercon- 
ducting transition T c fluctuations lead to the enhance- 
ment of the conductivity Q , while below T c they destroy 
the long- range order and lead to the finite resistance of 
the system. 

In the vicinity of T c for example in thin superconduct- 
ing wires with the diameter smaller than the coherence 
length thermal activation of the phase slips centers lo- 
cally destroys the superconductivity 0, S, @, [13] • Phase 
slip event is of the order of the coherence length £(T) 
where the amplitude of the order parameter vanishes at 
one point while the phase difference between the opposite 
sites of this point is ir. 

The observed resistance of the extremely thin super- 
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edge superconductivity still remains open. 

We will give the detailed analysis of the sample resis- 
tivity dependencies on temperature and magnetic field. 
The equation for Aslamazov- Larkin correction to the 
conductivity at fields higher than H c3 {T) will be also ob- 
tained. 

Let us consider a thin superconducting film with the 
magnetic field applied perpendicular to the surface of the 
film (see FigJTJ. We will study the case of the 2-type su- 
perconductor under the surface superconductivity con- 
dition, when the Ginzburg- Landau parameter k 1. 
Thus, we will not take into account the magnetic field 
modulations due to supercurrents. The size of the film 
is such that d <C £,(T) -C L, where d and L is the width 
and length of the film. The Ginzburg- Landau equation 
in dimensionless variables could be written in the form 



(zV + A) 2 f = *(1 - |*| 2 ) 



(1) 



where ^ is the complex order parameter, length is 
measured in units of the coherence length £(T) = 
(fnrD/8(T c (H) - T)) 1 / 2 , A- vector potential measured 



in units 
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D- diffusion coefficient. 
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conducting wires at temperatures T « T c 
argued to be caused by the phase slip events due to quan- 
tum tunneling of the order parameter [l3[ . 

In the present paper we will show that similar fluctua- 
tion effects appear in the edge superconductivity of thin 
superconducting film in the magnetic field perpendicular 
to the plane of the film. To the best of our knowledge 
the question of the fluctuations of the order parameter 
in the vicinity of the phase transition in thin films with 



Since the superconductivity in this regime exists only 
in the thin edge layer of the film we can treat each edge 
independently. Then let y axis to be applied along the 
corresponding edge, x axis directed to the bulk of the 
film. Magnetic filed is applied along the z-axis. Vector 
potential is chosen in the Landau gauge, A = (0, Hx, 0). 
The boundary condition for the order parameter at the 
(0, y) edge of the film is given as 

^L=o=0 (2) 

at the same time, the order parameter vanishes at the 
bulk of the film. We will search for the solution of the 
nonlinear Ginzburg- Landau equation |TJ in the form 

* (x,y)= 7 g(x)e ik v (3) 

where 7 is some constant, g(x) is a function subject to 
the condition J Q g 2 (x)dx — 1. Then the equation (JTJ 
can be reformulated as 



V(x) = \ lgo (x)\ 2 ,E = 



(4) 



2 



mate the solution for the g{x) by the function 
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FIG. 1: Thin film of thickness d in perpendicular magnetic 
field H , magnetic length Ih is the width of the superconduct- 
ing layer, L- the length of the layer along the edge of the 
film 



where go{x) is the eigenfunction corresponding to the 
lowest eigenvalue Eq of the equation 



ax 



1 + V{x) g(x) = Eg{x) (5) 



In order to solve equation ([5]) we can use the pertur- 
bation theory for 7 < 1. The eigenfunctions are then 
expressed through the parabolic cylinder function [l6| 

g n {x) cx D a (xV2H-k<f&), where a - E »- H+1 



2H 



The eigenvalue of the linear equation ([5]) corresponding 
to the ground state of the system is a function of fc and 
has a minimum where 



E a = lk-k (H)} 2 -e + 7 2 / 

Jo 



go(x)dx (6) 



where e = 1 — H/H c3 and k (H) ~ \Aff . First two terms 
in the equation |6j) can be obtained from the boundary 

^2H-k, 



condition -3-D, 



a \ Xy 6£1 - K y jjj \ x -0. 

The ground state eigenvalue condition min Eq = de- 
scribes the superconducting transition point. The choice 
of fc = ko(H) corresponds to the zero current case, fur- 
ther we will treat only this case. Finally, we obtain 
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dxgo(x) 



(7) 



We also provide the numerical solution for the nonlin- 
ear Ginzburg- Landau equation. Fig. [2] shows the zero 
current case k = ko(H) solution for different values of 
the magnetic field. Indeed, the order parameter is local- 
ized in the vicinity of the edge of the film at distance of 
the order of the magnetic length. The amplitude of the 
order parameter vanishes with increasing the magnetic 
field. The deviations from the solution obtained by the 
perturbation theory are small up to H ss H c2 . 

For practical purpose [l6| it is convenient to approxi- 



1/4 



(8) 



where b = y/l — 2/tt and H c3 = 1/6 ~ 1.66, that dif- 
fers from the exact one only by the 2%, while ko(H) = 
^H/bn. Using equations and §8§ we find the ap- 
proximate solution to the nonlinear Ginzburg -Landau 
equation 



*o(z,y) = (e\/2) 1/2 



(9) 



It is well known that phase slip events in ID supercon- 
ductin g w ires as well as vortices creation in superfiuid 
liquid T4j] at T < T c are due to thermal activation of 
the order parameter. The order parameter switches be- 
tween the metastable states of the superconductor chang- 
ing phase by 27r. The probability of such process is gov- 
erned by the Arrhenius law cx exp (-AF/T), where AF 
is the energy barrier separating these states. Solution for 
the phase slip event corresponds to the saddle point of 
the barrier. 

We have found the numerical solution for the phase 
slip center in the edge layer of the thin film by analyzing 
the time- dependent Ginzburg- Landau equation (for the 
review see [19J ) with the periodic boundary conditions 
on the order parameter. The solution for the amplitude 
of the order parameter is shown in Fig. [3] for different 
values of x at the magnetic field H = 1.6, k = 0.93, 
which corresponds to e = 0.1. In order to derive the 




FIG. 2: The amplitude of the order parameter |^o| for the 
H = 1, k = 0.73 (dashed line); H = 1.3, k = 0.85 (dash- dot); 
H = 1.6, k = 0.93 (solid line) 



approximate analytical solution for the phase slip cen- 
ter we will search for the order parameter in the form 
^i(x,y) = je lkaV Y /n >o C n(y)9n(x), where the summa- 
tion is over the set of eigenfunctions of the equation |5]| . 
In the zero-mode regime taking into account only n = 0, 
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we find 



cPC 
dy 



C 2 







(10) 



This equation is similar to the Ginzburg- Landau equa- 
tion for the one- dimensional wire at zero applied current 
case. Solving equation (fT0|) we obtain the expression for 
the order parameter corresponding to phase slip event in 
surface superconductivity 



&i(x,y) = # (z,y)tanh i 



(11) 



Taking into account solution (fTTj) we conclude accord- 
ing to [7[ with the expression for the resistivity of the 
thin film superconducting edge layer. 



R 



irh 2 Q 
2^T 



exp (-AF/T) 



(12) 



where 



= ^^fh tkA J dxdy (|* | 4 - |*i| 4 ) (13) 

is the saddle-point free energy barrier increment, f2 = 
(L/£(T))(e 3 / 2 /V GL ) (AF/T) 1/2 is the attempt frequency 
L is the length of the edge superconducting layer, tql — 
[ttTi i '8(T C 3(H) — T)] is the relaxation time and k is the 
GL parameter. 

Using equations © and QT]) at H < H c3 (T) we find 



AF — 



h -^fl[i H ^ H d]e^ (14) 



y/hc/eH c3 (T)- is 



where d- is the film thickness, £h c 
the magnetic length. 

This expression is simply the condensation energy of 
the superconductivity in volume £H o3 &Hd of thin film su- 
perconducting edge layer. 

Notice that the width of the edge superconducting 
layer £jj at T — > T C (H) is much smaller than the length 
of the normal part of the layer £h c3 /\^ caused by the 
thermal activation of the phase slip event, pointing the 
applicability of the thin wire approximation to the sur- 
face superconductivity of thin film. 

Phase slip events in ID superconducting wires at low 
temperatures are argued to be due to quantum tunnel- 
1, [13, EH- The resistivity of the wire is then 
exp (— 2S), where the exponent of the tunneling 
amplitude is S = A^G^ and A is a constant of the or- 
der of unity, is the conductance of the wire of length 
£(T) In our case of the edge superconducting layer 
it is the magnetic length £# c3 that governs the tunneling 
amplitude. The tunneling process is accompanied by the 
creation of the acoustic plasmons [13] which are responsi- 
ble for the interaction between the phase slip centers. As 
a result, this interaction suppresses the tunneling proba- 
bility and this effect is stronger as smaller the plasmons 
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FIG. 3: The amplitude of the order parameter for 
the case H — 1.6, k = 0.93; solid line corresponds to the 
solution on the boundary x = 0, dashed line corresponds to 
x = 1, x = 1.6, x — 2.3 



velocity, i.e. as stronger the Coulomb screening. The 
dissipation effects are the second factor that leads to de- 
creasing of the probability of the phase slip centers due 
to quantum tunneling [151 ]. 

In contrast to the case of isolated superconducting 
wire, the Coulomb interactions in the regime of edge su- 
perconductivity is screened by the charge of the normal 
part of the film. This immediately leads to decreasing of 
the acoustic plasmons velocity. The interaction with the 
normal part of the film results in the effective relaxation 
of the order parameter phase fluctuations. Consequently 
the resistivity of the surface superconducting layer of the 
thin film should be lower than the resistivity of the wire, 
taken under the same conditions. 

Let us consider the fluctuations of the order parameter 
in the case of surface superconductivity at magnetic fields 
H > H c3 (T). The corresponding Aslamazov- Larkin cor- 
rection to the conductivity at magnetic fields higher than 
the superconducting transition field H c % was studied in 
the paper [l8j] . It was shown that for the case of two- di- 
mensional surface superconducting layer this correction 
has the same temperature dependence as for the thin film 
at T > T c . 

However, in present work we focus on the extreme case 
of surface superconductivity, when the order parameter 
is concentrated in the quasi-one- dimensional layer of the 
thin film. 

Fluctuations of the order parameter at magnetic fields 
H > Hcz(T) result in additional correction to the con- 
ductance which according to [19( is 



G 



2m 2 



(15) 



where summation goes over the set v = n,k. The 
value of fluctuation of the order parameter ^S(x,y) = 



4 



Si/ 4>v9v{x)e v is written as 



(l^l 2 ) 



2m£ 2 Hc3 
bh 2 



T 



E n {k) 



(16) 



while t„ = tgl/ E n (k) is the characteristic decay time of 
the fluctuation. The lowest eigenvalue with n = gives 
the main contribution to the sum and coming from sum- 
mation to the integration over k we obtain the correction 
to the conductance of unit length of the layer 



G ~ 0.28 



e 2 H c3 (0) £ 

U H c3 (T) | e |3/2 



(17) 



Notice that the functional dependence of the correction 
G oc | e | 3/2 j s s i m ii ar to -t; ne case f one _ dimensional 
superconducting wire. 

It is seen also, with decreasing the temperature H c3 {T) 
increases and the value of the Aslamazov- Larkin correc- 
tion decreases. However, at the same time equation (|17p 
still valid in the interval of the magnetic fields H—H c3 (T) 
that increases with decreasing the temperature. 

For the numerical estimations we take typical values 
H c3 ~ IT then £h o3 ~ 25nm. For ultrathin film d ~ 
lOnm, k = 10, we estimate at T c3 ~ IK for the value 
AF/T ~ 10 3 e 3 / 2 . The probability of the phase slip event 
becomes negligibly small unless the parameter e = 1 — 
H/H c3 is of the order of 10" 2 . 

According [l6| the critical current destroying the sur- 
face superconductivity could be written as J c ~ jc(-H o3 d, 
where j c = — ^= is the critical current density of 

thin wire, we estimate J c ~ 20/iA 

To summarize, we have shown that the phase slip phe- 
nomenon reveals in edge superconducting layer of thin 
film in perpendicular magnetic field at H < H c3 (T). 
The corresponding resistance was calculated. The 
Aslamazov- Larkin correction to the edge superconduc- 
tivity of thin film at H > H c3 {T) have also been ob- 
tained. We conclude that such structures could be ap- 
plied as a new system for the study of the phase slip 
phenomenon in one- dimensional superconducting wires. 
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